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We investigate an evolutional scenario of the FRW universe with the Casimir energy scahng hke 
( — + The Casimir effect is used to explain the vacuum energy differences (its value measured 
from astrophysics is so small compared to value obtained from quantum field theory calculations). 
The dynamics of the FRW model is represented in terms of a two-dimensional dynamical system to 
show all evolutional paths of this model in the phase space for all admissible initial conditions. We 
find also an exact solution for non flat evolutional paths of Universe driven by the Casimir effect. The 
main difference between the FRW model with the Casimir force and the ACDM model is that their 
generic solutions are a set of evolutional paths with a bounce solution and an initial singularity, 
respectively. The evolutional scenario are tested by using the SNIa data, FRIIb radiogalaxies, 
baryon oscillation peak and CMB observation. We compare the power of explanation of the model 
considered and the ACDM model using the Bayesian information criterion and Bayesian factor. 
Our investigation of the information criteria of model selection showed the preference of the ACDM 
model over the model considered. However the presence of negative like the radiation term can 
remove a tension between the theoretical and observed primordial *He and D abundance. 

PACS numbers: 98.80.Jk, 04.20.-q 



I. INTRODUCTION 

The recent astronomical observations like SNIa P, i, WMAP i indicate that the current Universe is in an 
accelerating phase of its expansion. While there are many different explanations of this phenomenon fi*] the most 
straightforward one is that the Universe acceleration is due to the presence of a dark energy component comprising 
more than 70% of the energy of the Universe. Within this class of models the simplest candidate for dark energy is the 
phenomenological cosmological constant interpreted as vacuum energy Q. Unfortunately a simple estimation from 
quantum field theory gives that vacuum energy is larger than what we observe by a factor of the order 10^^". This 
discrepancy is called the cosmological constant problem and we are looking for some physical process which works to 
set vacuum energy precisely to zero. The Casimir force which is the manifestation of the quantum fluctuation can 
offer some mechanism to suppress the cosmological constant interpreted as vacuum energy to the value close to zero 
at the present epoch [1, 0j H Q- In this context the Casimir effect 10] seems to be relevant because it teaches us 
how vacuum energy contributes to the cosmological constant. Many other authors (see for example [Tlj ) argue that 
the Casimir effect is responsible for the vacuum energy differences between the value predicted by the quantum field 
theory and obtained from observations. Recently Antoniadis et al. [tI has suggested an important role of the Casimir 
effects in dark energy problem. 

Ishak pointed out the relevance of experiments at the interface of astrophysics and quantum field, focusing 



on the Casimir effects. Nontrivial topology of the Universe [12| was the motivations of these investigations in the 
cosmological context. Zeldovich and Starobinsky [l^ investigated the simple closed FRW models equipped in a 3- 
torus topology. The Casimir energy has been also studied in the context of compactification of extra-dimensions 
as an effective mechanism of dimensional reduction [l^ . ITsl flGj . It is worthy to mention that the Casimir type of 
contribution which arises from the tachyon condensation is also possible [l7j . 

For example in a toroidal model with the compactification scale L, one typically obtains a Casimir contribution 
p = (Tq) — —^a^ |18j] (see also ^9|]) with the scale factor a and the constant term a which depends on the nature 
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and the number of matter fields. The analogous result was obtained by Wreszinski who used the "cosmic box" idea 
of Harrison ^] to explain why the cosmological constant assumes an absurdally small value of energy density. For 
this aim local theory of the Casimir effect was applied to the Universe as a whole. The case of massless scalar field 
with no potential was also investigated in the cosmological context by Hardeiro and Sampaio [2]| . The authors 
study backreaction of the metric on quantum effects in terms of fluid with Casimir energy pcas = a/a^, pressure 
PCas — a/ia^ and the sign of a for a massless scalar as a function of coupling i^. The constant a changes the sign 
corresponds to a theory where the renormalized energy vanishes in the Einstein static universe. This critical value is 
Ccrit ^ 0.05391. 

The massless conformally coupled scalar field as well as the electromagnetic field and the massless Dirac field on the 
background of the Einstein static universe (or equivalently quantum effects calculated in the adiabatic approximations 
which can be justified by the fact that characteristic time scale of quantum process is smaller than a scale of time 
evolution) were considered in Refs. [l^jlHjH^I- In all cases the Casimir energy is of the form a/a^; when a is positive, 
than Casimir energy is an attractive force, while in an opposite case is repulsive [l3| . In the latter case Zeldovich 
and Starobinski [l3( suggested that the scalar field could drive inflation in a flat universe with a nontrivial-toroidal 
topology. It is interestiiig that a dynamical effect of the Casimir force is dynamically equivalent to effects of loop 
quantum gravity effects [25| . 

The FRW model with the Casimir type force contains a term which scales like negative radiation (— )(1 + z^. 
One should note that there are different interpretations of the presence of such a term: a cosmological model with 
global rotation, the Friedmann-Robertson- Walker (FRW) model in the Randall-Sundrum scenario with dark radiation, 
the FRW universe filled with a massless scalar field in a quantum regime or the FRW model in a semi-classical 
approximation of loop quantum gravity [2^ [27| , however in the present paper the FRW universe filled with a massless 
scalar field in a quantum regime (Casimir effect) is of our special interest. 

Our paper is organized as follows. In section II we provide some useful dynamical background for further presenta- 
tion. In section III we discuss observational constraints on the parameters of model under consideration. We complete 
the paper by Appendix in which we present an exact solution in the bouncing cosmology with the Casimir effect. 



II. DYNAMICS OF THE FRW MODEL WITH THE CASIMIR DARK RADIATION TERM 

In this section we investigate dynamical effects of the Casimir force which phenomenologically can be modelled by 
the Casimir energy pcas oc (— )a~^, where a is the scale factor. It can be shown that dynamical low temperature 
quantum effects of a scalar field can be modelled in terms of such perfect fluid [ij, [l^, [3l • It is useful to investigate 
the FRW dynamics with the cosmological constant and the Casimir force using dynamical system methods. The main 
advantage of this tool is possibility of investigations of all evolutional path for all admissible initial conditions. The 
main feature of cosmological models with the Casimir force is the presence of the bounce instead of the standard 
initial singularity. It is worth mentioning that it is not possible by using a geometrical test (basirig on null geodesies) 
to separate an individual Casimir component from all scaling in an analogous way, like (1 -I- z)'^ [26l. [27j. 

We consider the FRW model filled by a perfect fluid with energy p and pressure p. Then the evolution of the scale 
factor a can be described by the simple acceleration equation 

- = -^(p + 3p), (1) 
a 

where p and p are function only the scale factor. It is convenient to represent equation ([T]) in term of dimensionless 
variable x = where is the value of the scale factor at the present epoch in the form analogous to the Newtonian 
equation of motion, i.e. 

dV 

i = -^(2;), (2) 

where V{x) is the potential function and overdot here denotes the differentiation with respect to the rescaled cosmo- 
logical time t such that t ^ t: \Ho\dt = dr, Hq is the present value of the Hubble function. 

The Newtonian equation of motion ^ are equivalent to two-dimensional dynamical system of a Newtonian type 

dV 

x = y, y'^^'d^- 



Therefore, admissible critical points j/o = 0, = can be only saddles or centers. 

Let us consider the Universe with the cosmological constant {p\ = —pA, pa — A) filled by standard matter in which 
there are present the Casimir effect modelled as the phenomological fluid scaling like pcas oc (—)a^. The assumed 
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negative sign of energy pcas denote that the Casimir force is attractive. The potential is 

V{x) = -^{flAfiX^ - Ofc,o + i^m,OX^^ + flci^sfiX-^), (4) 

where r is new time parameter proportional to the original cosmological time t such that t t: (HQ)dt = dr; 
X = a/ao is a dimensionless scale factor expressed in the units of its present value gq, f2i^o are the density parameters 
for matter (i m), the Casimir contribution {i = Cas), the cosmological constant (i — A) and the curvature k (the 
curvature index fc = 0, ±1). 

Equation (j4]) admits the first integral 

U_ + v{x)=E^O (5) 

which has a simple interpretation of conservation energy. The universe is presently {x = 1) accelerating if Vix) is a 
decreasing function of its argument x: (^)^_]^ < 0. 

From equation ([5]) we obtain at the present epoch the constraint V (x = \) = —\ oy 

I] a- 0-1- (6) 

i 

From equation Q we obtain that negative values of a as well as the positive cosmological constant give rise to 
acceleration of the Universe. Because of constraint ([5]) the motion of the system takes place in a domain Vq = 
x: V{x) < 0) of the configuration space. 

If we shift the curvature term —■^flk,o on the right hand side of ([5]) then system ^ can be considered on the 
constant energy level E = — irj^ o parametrized by the curvature constant. The phase space {x,x') of all trajectories 
for all admissible initial conditions is divided by trajectory of the flat models flk,o = 0. Closed models {ilk,o < 0) are 
situated inside this characteristic curve and open models {^k.o > 0) outside. 

The phase portraits together with the potential function ([4]) are shown in Fig. [l]and Fig. [2] for a < and a > 0, 
respectively. From the phase portrait in Fig. [2] (for a < 0, pcas = aa~^) we can observe that the initial singularity 
characteristic for the case of a > is replaced by a bounce. To have a bounce, there must be some time at which 
the size of the universe (or the scale factor) assumes a minimum: a{to) = 0. Moreover d{to) > is required. The 
sufficient condition for the origin of the FRW universe from a bounce was formulated by Molina-Paris and Visser ^2§\ . 
From the point of view of the structural stability notion there is the main difference between phase portraits in Fig. [1] 
and Fig. [21 The phase portrait for the ACDM model is structurally stable. This means that small perturbation of 
right-hand sides does not disturb the phase portrait. In contrast to Fig. [T] the phase portrait in Fig. [21 is structurally 
unstable due to presence of centre on x-axis. Following the Peixoto theorem the ACDM dynamical system is generic in 
the class of all planar dynamical systems and bouncing cosmological models are exceptional because do not form such 
open dense subsets. The type of a critical point is determined by eigenvalues of the linearization matrix calculated at 

the critical point. This equation assumes the form -I- — 0. Therefore the bounce representing by a centre 

X—Xq 

is present on the phase portrait if Vxxixo) > 0, where Vx{xo) = 0, xq is the critical point geometrically it means that 
V{x) function has minimum in xq. 

Note that all models on the phase portraits behaves reflectional symmetry x — > {—)x, x' — > {—)x'. During the 
bounce the SEC (strong energy condition) violation is necessary. 

For present acceleration of the universe it is required 

— ^m,0 — 2r2cas,0 + f^A,0 > 0. (7) 

Therefore if A = the Casimir force must be sufficiently large 

(|f^Cas,0|+f^A,0)>%^ (8) 

The bouncing cosmology is following of Barrow's reincarnation of ancient fascination of cyclicity and realization 
myth of the "eternal return" [29,]. While the classical bounces are generated by fields endowed with negative energy, 
the quantum gravity effects in semi-classical approximation invoked to justify the bounce [s^l generated by some 
quantum corrections terms. Also possibility that a contracting braneworld model experiences a bounce instead ever 
reaching a singularity is addressed by Shtanov and Sahni [sil . [3^ . Issl . IsT 



Recently bouncing cosmology has been analyzed in details in 
a bounce is close to unity for models of f{R) nonlinear gravity 



It was also demonstrated that the probability of 

[13, Hi. 
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FIG. 1: Phase portrait for the FRW model with a > (usual radiation) and the cosmological constant A (positive). There are 
three different types of evolution oscillating, loitering and bouncing. The trajectory fc = of the flat model separates all the 
evolutional paths of closed and open models. 



III. TOWARD TESTING BOUNCING COSMOLOGY CAUSED BY ( )(1 + TERM 

It is very interesting to test cosmological models against observations. One of the most popular test used the 
SNIa data. This test is based on the luminosity distance of the supernovae la as a function of redshift The 
observations of the type la distant supernova suggest that the present Universe is accelerating P, 0, [H, HO] . Every 
year new SNIa enlarge the available data by more distant objects and lower systematics errors. Our work is based 
on two samples: Riess et al. [s^ "Gold" sample of 157 SNIa and Astier et al. [4l| sample of supernovae, based 
on 71 high redshifted SNIa discovered during the first year of the 5-year Supernovae Legacy Survey. We used this 
sample because we would like to easily cornpare the result obtained in present paper with that obtained, with using 
of different method, in our previous paper [26|. 

One should note that for the distant SNIa, we directly observe their apparent magnitude m and redshift z. The 
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FIG. 2: Phase portrait for the FRW model with the repulsive Casimir force (a < 0). All models are undergoing a bounce 
instead of an initial singularity replaced by the critical point {x,x) — (0, co) like in Fig. 1. In the phase portrait there are 
present two types of static critical points center and saddle point like in Fig. 1. Note that bounce appears to be possible in 
open flat and closed cosmologies. The phase portrait is structurally unstable in contrast to the ACDM model represented in 
Fig. 1. 

absolute magnitude Ai of the supernovae is related to its absolute luminosity L. We have the following relation 
between distance modulus /x, the luminosity distance (11, the observed magnitude m and the absolute magnitude M: 

fi = m- M = 5 logio + 25 = 5 logio Dl+M (9) 

where Dl ~ HqcIl and A4 = —5 logj^p Hq + 25. The luminosity distance of a supernova is the function of cosmological 
parameters and redshift 
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FIG. 3: The 68.3% and 95.4% confidence levels obtained from combined analysis SN+RG) on the (r2m,o,fiA,o) plane. 



where 



H 
i5n 



(11) 



^kM — and ^-"(2:) = (sinh(a;), sin(a;)) for fc < 0, fc = 0, A: > 0, respectively. We assumed f2r,o = ^^7,0 + ^^1^,0 — 

2.48/1-2 X 10^5 _^ ij/1-2 X 10-5 ~ 0.0001 m. 

Substituting (fTU)) back into equations ^ provides us with an effective tool (the Hubble diagram) to test cosmological 
models and to constrain their parameters. Assuming that supernovae measurements come with uncorrelated Gaussian 
errors, one can determine the likelihood function C from chi-square statistic C oc exp(— x^/S), where 



E 



(m: 



thcor ,,obs^2 



obsVI 
H-i ) 



(12) 



The probability density function (PDF) of cosmological parameters [l| can be derived from Bayes' theorem. Therefore, 
one can estimate model parameters by using a minimization procedure. It is based on the likelihood function as well 
as on the best fit method minimizing x^. Constraints for the cosmological parameters, can be obtain by minimizing 
the following likelihood function C oc exp(— y2/2). 



Daly and Djorgovski [43| (see also [26 



|45[) suggested to include in the analysis not only supernovae but also 
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FIG. 4: The 68.3% and 95.4% confidence levels obtained from combined analysis SN+RG+BOP) on the (f^m.o, i^A.o) plane. 



radio galaxies. In such a case, it is useful to use the coordinate distance defined as [46 

HodL{z) 
oil + z) 

In such a case we can determinate likelihood function C oc cxp(— where 



E 



obs 



th 



(13) 



(14) 



To analyse SNIa data and FRIIb data in the unique way it is useful analises for both sets of data coordinate distance 
y{z). Daly and Djorgovski [i^] compiled a sample comprising the data on y{z) for 157 SNIa in the Riess et al. Isol 
Gold dataset and 20 FRIIb radio galaxies. In our data sets we also include 115 SNIa compiled by Astier et al. [4ll |. 
In our previous paper [261] there are more details on the computation of y{z) and (Ji{yi) for these samples. 

In Fig. Owe present the 68.3% and 95.4% confidence levels obtained from combined analysis SNIa and FRIIb on 
the plane (rim^O: ^^A.o)- 

We also include to our analysis the baryon oscillation peaks (BOP) detected in the Sloan Digital Sky Survey (SDSS) 
luminous red galaxies [481 . They found that the value A is 



A: 



■\/f^m,0 / 1 



k,0\ 



dz 



0.469 ±0.017 



(15) 
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The 68.3% and 95.4% confidence levels obtained from combined analysis SN+RG+BOP+CMB) on the (f2m,o, ^^a,o) 



where E{z) = H{z)/Hq and zi — 0.35. The quoted uncertainty corresponds to one standard deviation, where a 
Gaussian probability distribution has been assumed. These constraints could also be used for fitting cosmological 

" A'""'' -0.469 ' 



parameters [41l |49| . In such a case we determinate the likelihood function C cx exp 



0.017 



/2 



Another constraint which we also include in our analysis is the so called the (CMBR) "shift parameter" [5C 




dz 

W) 



= 1.716 ±0.062 



which leads to the likelihood function C cx exp 

In Fig. |4]and Fig. [5] we present confidence levels obtained from the analysis of SNIa, FRIIb and BOP or BOP and 
CMB shift, respectively. One should note that when we are interested in constraining a particular mod el p arameter, 
the likelihood function marginalized over the remaining parameters of the model should be considered [51| . 

From our combined analysis (SN+RG+SDSS+CMBR) we obtain as the best fit a flat (or nearly flat universe) 
with Om,o — 0-3, and r^A.o — 0-7. For the dark radiation term, we obtain the stringent bound r^R^o = ^^r.o ~ 
l^^dr.ol > —0.00025 at the 95% confidence level where fidr,o is negative. It leads for the limit on dark radiation 
l^^dr.ol < 0.00035. This results mean that null value of dark radiation term is preferred (fidr.o — 0), however the small 
negative contribution of dark radiation is also available. Our results shows that in the present epoch contribution of 
the dark radiation, if it exist, is small and gives only small corrections to the ACDM model in the low redshift. 
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Our result is in agreement with the result of our previous paper '2u\ where in the combined analysis based in a 
pseudo-x^ merit function (sij : we obtain as the best fit a flat universe with flni,o — 0-3i ^dr,o = and ^a,o — For 
the dark radiation term, we obtain the stringent bound |ildr,o| < 0.00035 at the 95% confidence level (|r2dr,o| < 0.00026 
at the 68.3% confidence level). 

Please also note that if SIr^o = ^i.o ^ l^^dr.ol < 0, then wc obtain a bouncing scenario [H, [H, [s^l instead of a big 
bang. For ilm.o — 0.3, ridr,o = —0.00035 and h = 0.65 bounces {H^ = 0) appear for z ~ 1200. In this case, the BBN 
epoch never occurs and all BBN predictions would be lost. 

We use the Akaike information criterion (AIC) [H^, the Bayesian information criterion (BIG) ^5^] as well as the 
Bayes factor to select model parameters providing the preferred fit to data. The information criteria put a threshold 
which must be exceeded in order to assert an additional parameter to be important in explanation of a phenomenon. 
The discussion how high this threshold should be caused the appearing of many different criteria. The Akaike and 
Bayesian information criteria (AIC and BIC) (for review see |55|) are most popular and used in everyday statistical 
practices. The usefulness of using the information criteria of model selection was recently demonstrated by Liddle 
[56l | and Parkinson et al. [53| • The problem of classification of the cosmological models on the light of the information 
criteria on the base of the astronomical data was discussed in our previous papers [H, [6ll, [H, [fi^l ■ 

The AIC is defined in the following way 

AIC = -21n/: + 2d (16) 

where C is the maximum likelihood and d is a number of the free model parameters. The best model with a parameter 
set providing the preferred fit to the data is that minimizes the AIC. The BIC introduced by Schwarz ^] is defined 
as 

BIC = -21n£ + dln7V (17) 

where N is the number of data points used in the fit. Comparing these criteria, one should note that the AIC tends 
to favor models with large number of parameters unlike the BIC, because the BIC penalizes additional parameters 
more strongly. Of course only the relative value between the BIC of different models has statistical significance. The 
difference of 2 is treated as a positive evidence (and 6 as a strong evidence) against the model with the larger value 
of the BIC [U, [6^1 . If we do not find any positive evidence from the information criteria the models are treated as a 
identical and eventually additional parameters are treated as not significant. 

In the Bayesian framework quality of the models can be compared with help of evidence [63, [6^ . We can define 
the a'posterior odds for two models Mi and Mj - the Bayes factor Bij [g^l . If we do not favor any model it reduces 
to the evidence ratio. Schwarz (s^ showed that for observations coming from a linear exponential family distribution 
the asymptotic approximation — > oo the logarithm of evidence is given by 

ln.B = ln/:-^lniV + 0(l). (18) 

It is easy to show that in this case we have the simple relation between the Bayes factor and the BIC 

21nBy = -(BICi - BICj) (19) 

If Bij is greater than 3 it is called positive evidence in favor of Mi model, while Bij > 20 give strong and Bij > 150 
very strong evidence in favor of Mi model (6^ . 

Our results are presented in Table HI In our investigations we consider the likelihood function C oc exp(— x^/2) 
where a pseudo-x^ merit function is used [l^, [13, [5l| : 

2_2 |2 |2 _ 

X — XSN+RG + XSDSS + XCMBR — 



^\ <y.) J \ 0.017 ) ^\ 0.062 ) ' ^'"^ 

where A^°'^ and denote the values of A and R obtained for a particular set of the model parameter. 

Our investigations of the information criteria show that the bouncing A cold dark matter model (BACDM) with 
dark energy does not increase the fit significantly. It confirms our conclusion that the dark energy term, if it exists, 
is small in the present epoch. The Bayes factor also favors the ACDM model over the BACDM model with dark 
energy. Its supports results obtained with help of the AIC and BIC. One should note that because of a non-Gaussian 
distribution of the a'posteriori PDF function for VIr.q the above results obtained with help of the Bayes factor should 
be treated with caution as only an additional support for results obtained with help of the AIC and BIC. 
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IV. CONCLUSION 



In this paper we analyzed the observational constraints on the (l + z)'*-type contribution in the Friedmann equation. 
The analysis of SNIa data as well as both SNIa and FRIIb radio galaxies with the constraints coming from baryon 
oscillation peaks and CMBR "shift parameter" shows that influence of a negative term —(1 + z)^ is very weak in the 
present epoch of the Universe. 

The main aim of this paper is investigation of evolutionary paths of the FRW cosmological models with the 
cosmological constant A and the Casimir energy scaling like (— )(1 + z)^- We demonstrate that a bounce is a generic 
feature of this model. In this scenario an initial singularity is replaced by a bounce. We characterize a full class of all 
evolutionary paths by using dynamic system methods and we find exact solution. We pointed out the main difference 
between the ACDM model and the BACDM model from the view of point structural stability. The two-dimensional 
case of dynamical system is distinguished by the fact that the Peixoto theorem gave a complete characterization of 
the structurally stable systems on any compact space, assert that they form open and dense subsets in the space 
of all dynamical systems on the plane. We conclude that while the ACDM models are structurally stable, therefore 
typical, the BACDM models are structurally unstable, therefore exceptional. Of course structurally unstable models 
may model a realistic physical situation however they require a fine tuning of the model itself (not fine tuning of 
initial conditions). 

We show that there are several interpretations of the (l-|-z)'*-type contribution and we discussed different proposals 
for the presence of such a term. Unfortunately, it is not possible, with present kinematic astronomical tests, to 
determine the energy densities of individual components scales like radiation. However we show that some stringent 
bounds on the value of this total contribution can be given. The combined analysis of SNIa data and FRIIb radio 
galaxies, baryon oscillation peaks and CMBR "shift parameter" give rise to the concordance universe model which 
is almost flat with Qrn,o — 0.3. From the above-mentioned combined analysis, we obtain an constraint for the 
negative term which scales like total radiation r^R^o > —0.00025 which leads to bounds on the dark radiation term 
^^dr,o > —0.00035 at the 95% confidence level. This is in agreement with result of our previous paper |26] obtained 
with help of a pseudo-x^ merit function. 

The investigations of the information criteria show the BACDM model with dark energy does not increase the fit 
significantly. It confirms our conclusion that the dark energy term, if it exists, is small in the present epoch. 

In this paper we have been especially studied the advance of an initial singularity using back reaction gravity 
quantum effect at low temperatures (the Casimir effect). The Casimir force arising from the quantum effect of 
massless scalar field give rise to a (— )(1 -I- 2;)^ correction whose effect depends upon the geometry and nontrivial 
topology of the space. Typically this type of correction is thought to be important at the late time of evolution of 
the universe. We have shown that the Casimir effect could generically remove the initial singularity which would be 
replaced by the bounce. 

Our analysis of back reaction on quantum effect clearly reflects an important role played by vacuum energy. The 
basic problem which helps us to understand the dynamics of accelerating universe is how the cosmological constant 
(geometrical) contributes to the vacuum energy. 
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The Friedmann equation we are dealing with can be reduced so that the right-hand side is a 4th degree polynomial. 
Such equations have been completely analyzed in 68] , so we merely point out which cases of the general classification 
(by which we will mean the cited paper's content) apply and what additional features arise here. 

The equation reads: 
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APPENDIX A: BASIC MODELS WITH THE CASIMIR FORCE. 




(Al) 



Analogously to the f2fl,o in section III we introduce Qji = fi^.o — smd change the variables to: 



y — X, du = xdr, 
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TABLE I: The values of AIC and BIC and Bayes factor S12 for ACDM model and bouncing cosmology model (with dark 
radiation). The upper section of the table represents the constraint fi^.o = (flat model). 





ACDM 


BACDM 




sample 


AIC BIC 


AIC BIC 


B12 


SN 


299.5 303.1 


301.5 308.7 


16.44 


SN+RG 


322.4 326.1 


324.4 331.8 


17.29 


SN+RG+SDSS 


324.4 328.1 


326.4 333.8 


17.29 


SN+RG+SDSS+CMBR 


324.5 328.2 


326.5 333.9 


17.29 


SN 


300.0 307.2 


302.0 312.8 


16.44 


SN+RG 


323.5 330.9 


325.5 336.5 


17.29 


SN+RG+SDSS 


325.1 332.5 


327.1 338.1 


17.29 


SN+RG+SDSS+CMBR 


326.5 333.9 


328.5 339.5 


17.29 



the main equation becomes 



^^A.oa;'* - ^k,Ox'^ + ^ImfiX + fifl = W{x). 



(A2) 



The only new solutions appear, when one considers iluj big enough, for flu to be negative. Otherwise, the solutions 
is equivalent to that of a model with radiation only. We note that: 



applies to all models, restricting the solutions. 

The most general form of the solution, satisfying a;(0) = 1, reads 



(A3) 



1 



p'{t) + Oa,o 

2[p(r) - ^(60a,o - ^^fc,o)]2 - iJ^A,o 



+ 



+ 



with the invariants 



52 — J2^k,0 + ^A,O^I 



93 — 2T6^fc,0 ~ ]^^A,0^^m,0 ^ l^A.O^kfl^R- 



(A4) 



(A5) 



and using it, we can also obtain an explicit formula connecting the conformal and cosmological times. 

First, we need to investigate the singular points of our solution. As x(t) is an elliptic function of order two, there 
are exactly two poles (provided that Oa,o 7^ 0). To specify them exactly, we can use the p function, together with its 
derivative 



Pin) 

P{t2) 



Relation du = xdr then becomes 



p'{ri) = +\^/nZo 
P'(r2) = -4 



A,0, 



1 



A,0 



= [C(Ti)-C(T2)]^r + 



:ln 



Jx/^^A.O + ^^A,0- 

a{T - n)a{T2) 



a{T - T2)<j{ti] 



(A6) 



(A7) 



where the constant of integration was choosen so that w(r = 0) = 0, and C, and a are the appropriate Weierstrass 
functions. 

This formula remains true in the degenerate cases as well, with C, and a also simplified. In the case of trigonomet- 
ric/hyperbolic solutions we have 



C(r) 



.gr + ^/Qg coi{J \g t) 



[2 1 2 
cr(r) = 4 / — e4^'^ sin(-^ 

V 35 



i^r), 



(A8) 
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FIG. 6: A generic oscillating solution, which in the considered case is never singular. Scaling from t to u only results in 
stretching the time axis. 

with \/3g = ±v^, and for initial conditions chosen as for the general solution. 

The two following sections consist of the case by case classification analogous to that of [68*] with the quantities 
ai,W-,x+f-fo defined therein. 

1. !^A,o > 
1.1. Case ae / 

Analyzing the behavior of ag, with the possible values of and f^A^Oi we obtain the first important restriction: 

nk,o sC ^ (76 < 0. (A9) 

For positive curvature, any sign is possible. 

As follows from the general considerations, we have two possibilities: Cg > 0, or ct@ < 0. However, the former is 

limited to the case of four real roots. This happens because we have: W- — {2Qk,o ~ y^il^, q + 12f7i^i7A,o)/3, and 

by (|A9[) rife.o > while Sl/friA.o < 0, causing w_ to be always positive, and making the resolvent cubic admit three 
solutions. 



1.1.1 Four complex roots 



This case is impossible, because JIa.o > would force w to be positive everywhere. 
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FIG. 7: A typical bounce in finite time r, corresponding to u changing from — oo to oo. 



1.1.2. Four real roots 



There are two possibilities here. Either ei < < 62 < 63 < 64, or ei < 62 < 63 < < 64. The former admits both 
the solutions xq and a;+, the latter only However, none of them passes through zero, allowing for the particularly 
interesting non-singular, periodic behavior of a;o. A typical oscillatory solution is presented in Fig.[6l and the bounce 
in Fig. [7l Although the latter has a finite time of evolution in r, this interval is stretched into the whole real line 
when in the variable u. This happens because the right hand side in (jA7p is singular. 



1.2. (T6 = — one double root 



As a consequence of (|A9p . this, and obviously all the following cases, require positive curvature. Also, as in the 
general classification, we must have: 



(AlO) 



As fjg — 0, is effectively an equation of degree two in f^A^o it instantly provides new restrictions. First, for it to have 
real solutions in flA,o, we must have 



and if at least one of the solutions is to be positive, we also need 

lefiflrife.o , 



(All) 



finally yielding: 



"m,0 



fc,0 



02 

"m,0 



3(3 + V3);-3(3- V3) 



-|;-3(3- V3) 



(A12) 
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FIG. 8: The quasi-static solution, which stays almost at the same value of the scale factor for most of the time (both r and u), 
and undergoes only one short dip. 



This also allows us to see that the sign oia^, which is the same as that of the expression n\^Q — 2k^ / {9flf-^^ Q+Silfc^o^^i?) 
in this case, can be both plus and minus. The biggest restriction comes from the inequality (IA3p . which makes case 
1.2.1 impossible. As for 1.2.2, all three subcases might happen. If the double root is the smallest, then we only have 
the solution (bounce). If it lies between the other roots, we have both the static, stable solution xq, and a;+ 
(bounce). Lastly, if the double root is the biggest, we have both asymptotic branches x+ and X-i, together with the 
unstable, static xq. 

The bouncing solutions are essentially the same as the previous ones (Fig. [7]) , also with time stretching. As for the 
asymptotic solutions, we either have a quasi-static evolution x^ with one minimum, depicted in Fig.[51 or a monotonic 
expansion/contraction presented in Fig. [5) The former has an infinite evolution time in both r and u, and the 
latter reaches infinite x for finite r but infinite u. 



1.3. (Te = (Ts = 0, (T4 7^ — two double roots 

Further simplifications require 

o "fc,o 
"A.o 



i2n 



R 



^m,0 = yJ-fnkfi^R. (A13) 

However, such a value of Qa.o makes = 0, making this case impossible, and bringing us to the next one. 



1.4. (76 = (75 = (74 = 0, (7a 7^ — One triple root 

The so far excluded possibility of 12riflSlA^o + o ~ 0' holds here, and substituting (|A13p into 0-3, we see that 
further simplification is not possible, as (T3 is proportional to ^ 0. Following the general classification 1.4, we 
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FIG. 9: Only the increasing solution in this figure is physical, it tends to the static solution at minus infinity (both r and u) 
and the expansion to infinity takes infinite time (u). 



can see that the multiple root must be negative, because of (jA3p the simple root must be positive, and thus only the 
bounce solution can be physical. Again, Fig. [7] applies. 



2. Qa.o < 



2.1. (76 7^ — simple roots 



As before we obtain a restriction on fTg 



^k,0 < ^ (76 > 0, 



(A14) 



and similarly to before, for non-negative curvature any sign is possible. Analyzing the subcases of erg > 0, we can see 
that one of the restrictions becomes an identity now. Namely: 



k.O ■ 



(A15) 



is always true. Furthermore, for A; ^ 0, we can only have w_ < 0, but for A; > both signs are possible. In 
consequence, we can have the subcase 2.1.2, which is all the more interesting thanks to (jA3[) . which ensures, that 
there are oscillations with no singularity; and 2.1.3 with two possible regions of oscillation, also without singularity. 
They are qualitatively the same as the > oscillations depicted in Fig. [51 



2.2. (76 = — one double root 



Because of (IA14p . from now on we must have non-negative curvature. Obviously, (|A11[) must still hold, but the 
intervals of flj^ are now "reversed" : 



e (-cx);-3(3 + V3) U -3(3-V3);0), 
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finally giving: 




e -3(3- \/3);0) . 



(A16) 



Looking at (T5 as a function of Qa.o with the above restrictions, it is straightforward to check that it can only be 
negative, thus reducing this case to 2.2.1 only, where only a static, stable solution exists. 



As CTs = has only positive roots with all the current restrictions on the parameters, this and all the further cases 
are impossible. 



[1] A. G. Riess et al. (Supernova Search Team), Astron. J. 116, 1009 (1998), astro-ph/9805201. 

[2] S. Porlmuttor ot al. (Supernova Cosmology Project), Astrophys. J. 517, 565 (1999), astro-ph/9812133. 

[3] D. N. Spergel et al. (WMAP), Astrophys. J. Suppl. 148, 175 (2003), astro-ph/0302209. 

[4] E. J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. Phys. D15, 1753 (2006), hep-th/0603057. 

[5] T. Padmanabhan, Phys. Kept. 380, 235 (2003), hep-th/0212290. 

[6] M. Ishak (2005), astro-ph/0504416. 

[7] 1. Antoniadis, P. O. Mazur, and E. Mottola, New J. Phys. 9, 11 (2007), gr-qc/0612068. 
[8] M. Bordag, U. Mohideen, and V. M. Mostepanenko, Phys. Rept. 353, 1 (2001), quant-ph/0106045. 
[9] S. A. EUingscn and I. Brcvik, J. Phys. A40, 3643 (2007), quant-ph/0611030. 
[10] H. B. G. Casimir, Kon. Ned. Akad. Wetensch. Proc. 51, 793 (1948). 

[11] K. A. Milton, R. Kantowski, C. Kao, and Y. Wang, Mod. Phys. Lett. A16, 2281 (2001), hep-ph/0105250. 

[12] M. Lachieze-Rey and J.-P. Luminet, Phys. Rept. 254, 135 (1995), gr-qc/9605010. 

[13] Y. B. Zeldovich and A. A. Starobinsky, Sov. Astron. Lett. 10, 135 (1984). 

[14] M. Szydlowski, J. Szczesny, and M. Biesiada, Class. Quantum Grav. 4, 1731 (1987). 

[15] M. Szydlowski and J. Szczesny, Phys. Rev. D38, 3625 (1988). 

[16] M. Szydlowski, J. Szczesny, and T. Stawicki, Class. Quant. Grav. 5, 1097 (1988). 

[17] B. Mclnnes (2006), hep-th/0607074. 

[18] V. M. Mostepanenko and N. N. Trunov, The Casimir Effect and its Applications (Clarendon Press, Oxford, 1997). 

[19] D. Muller (2004), gr-qc/0403086. 

[20] W. F. Wreszinski (2006), quant-ph/0603026. 

[21] C. A. R. Hcrdeiro and M. Sampaio, Class. Quant. Grav. 23, 473 (2006), hep-th/0510052. 
[22] J. S. Dowker and R. Critchley, Phys. Rev. D15, 1484 (1977). 
[23] E. Strccruwitz, Phys. Rev. Dll, 3378 (1975). 

[24] F. W. J. Olver, Asymptoties and Special Functions (Academic Press, New York, 1974). 

[25] D. J. Mulryne, R. Tavakol, J. E. Lidsey, and G. F. R. Ellis, Phys. Rev. D71, 123512 (2005), astro-ph/0502589. 

[26] W. Godlowski and M. Szydlowski, Phys. Lett. B642, 13 (2006), astro-ph/0606731. 
[27] W. Godlowski, M. Szydlowski, and Z.-H. Zhu (2007), astro-ph/0702237. 
[28] C. Molina-Paris and M. Visser, Phys. Lett. B455, 90 (1999), gr-qc/9810023. 

[29] J. D. Barrow, D. Kimberly, and J. Magueijo, Class. Quant. Grav. 21, 4289 (2004), astro-ph/0406369. 
[30] P. Singh and A. Toporensky, Phys. Rev. D69, 104008 (2004), gr-qc/0312110. 
[31] Y. V. Shtanov (2000), hep-th/0005193. 

[32] Y. Shtanov and V. Sahni, Class. Quantum Grav. 19, LlOl (2002), gr-qc/0204040. 
[33] V. Sahni and Y. Shtanov, JCAP 0311, 014 (2003), astro-ph/0202346. 
[34] Y. Shtanov and V. Sahni, Phys. Lett. B557, 1 (2003), gr-qc/0208047. 

[35] M. Szydlowski, W. Godlowski, A. Krawiec, and J. Golbiak, Phys. Rev. D72, 063504 (2005), astro-ph/0504464. 

[36] A. Borowiec, W. Godlowski, and M. Szydlowski, Phys. Rev. D74, 043502 (2006), astro-ph/0602526. 

[37] A. Borowiec, W. Godlowski, and M. Szydlowski, Int. J. Geom. Meth. Mod. Phys. 4, 183 (2007), astro-ph/0607639. 

[38] S. Carloni, P. K. S. Dunsby, and D. M. Solomons, Class. Quant. Grav. 23, 1913 (2006), gr-qc/0510130. 
[39] A. G. Riess ct al. (Supernova Search Team), Astrophys. J. 607, 665 (2004), astro-ph/0402512. 
[40] A. G. Riess ct al. (2006), astro-ph/0611572. 

[41] P. Astier et al., Astron. Astrophys. 447, 31 (2006), astro-ph/0510447. 

[42] R. G. Vishwakarma and P. Singh, Class. Quant. Grav. 20, 2033 (2003), astro-ph/0211285. 
[43] R. A. Daly and S. G. Djorgovski, Astrophys. J. 597, 9 (2003), astro-ph/0305197. 
[44] Z.-H. Zhu, M.-K. Fujimoto, and X.-T. He, Astrophys. J. 603, 365 (2004), astro-ph/0403228. 
[45] D. Puetzfeld, M. Pohl, and Z.-H. Zhu, Astrophys. J. 619, 657 (2005), astro-ph/0407204. 
[46] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 



2.3. ae = (75 = 0, (74 ^ — two double roots 



17 



R. A. Daly and S. G. Djorgovski, Astrophys. J. 612, 652 (2004), astro-ph/0403664. 
D. J. Eisenstcin ct al., Astrophys. J. 633, 560 (2005), astro-ph/0501171. 
M. Fairbairn and A. Goobar, Phys. Lett. B642, 432 (2006), astro-ph/0511029. 
Y. Wang and M. Tegmark, Phys. Rev. Lett. 92, 241302 (2004), astro-ph/0403292. 

V. F. Cardone, C. Tortora, A. Troisi, and S. Capozziello, Phys. Rev. D73, 043508 (2006), astro-ph/0511528. 
B. K. Tippett and K. Lake (2004), gr-qc/0409088. 
H. Akaike, IEEE Trans. Auto. Control 19, 716 (1974). 

G. Schwarz, Annals of Statistics 6, 461 (1978). 

K. P. Burnham and D. R. Anderson, Model Selection and Multimodel Inference: A Practical Information - Theoretical 
Approach (Springer- Vcrlag, New York, 2002), 2nd ed. 

A. R. Liddle, Mon. Not. Roy. Astron. Soc. 351, L49 (2004), astro-ph/0401198. 

D. Parkinson, S. Tsujikawa, B. A. Bassett, and L. Amendola, Phys. Rev. D71, 063524 (2005), astro-ph/0409071. 
W. Godlowski and M. Szydlowski, Phys. Lett. B623, 10 (2005), astro-ph/0507322. 
M. Szydlowski and W. Godlowski, Phys. Lett. B633, 427 (2006), astro-ph/0509415. 
M. Szydlowski and W. Godlowski, Phys. Lett. B639, 5 (2006), astro-ph/0511259. 
M. Szydlowski and A. Kurek (2006), gr-qc/0608098. 

M. Szydlowski, A. Kurek, and A. Krawiec, Phys. Lett. B642, 171 (2006), astro-ph/0604327. 
A. Kurek and M. Szydlowski (2007), astro-ph/0702484. 

H. Jeffreys, Theory of Probability (Oxford University Press, Oxford, 1961), 3rd ed. 

S. Mukherjee, E. D. Feigelson, G. J. Babu, F. Murtagh, C. Praley, and A. Raftery, Astrophys. J. 508, 314 (1998), astro- 

ph/9802085. 

P. Mukherjee, D. Parkinson, P. S. Corasaniti, A. R. Liddle, and M. Kunz, Mon. Not. Roy. Astron. Soc. 369, 1725 (2006), 
astro-ph/0512484. 

R. E. Kass and A. E. Raftery, J. Amer. Stat. Assoc. 90, 773 (1995). 

M. P. Dabrowski and T. Stachowiak, Annals Phys. 321, 771 (2006), hep-th/0411199. 



